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The  goal  of  this  work  is  to  develop  a  numerical  model  for  the  vibration  of  hybrid  electric  vehicle  (HEV) 
battery  packs  to  enable  probabilistic  forced  response  simulations  for  the  effects  of  variations.  There  are 
two  important  types  of  variations  that  affect  their  structural  response  significantly:  the  prestress  that  is 
applied  when  joining  the  cells  within  a  pack;  and  the  small,  random  structural  property  discrepancies 
among  the  cells  of  a  battery  pack.  The  main  contributions  of  this  work  are  summarized  as  follows.  In 
order  to  account  for  these  two  important  variations,  a  new  parametric  reduced  order  model  (PROM) 
formulation  is  derived  by  employing  three  key  observations:  (1)  the  stiffness  matrix  can  be  parame¬ 
terized  for  different  levels  of  prestress,  (2)  the  mode  shapes  of  a  battery  pack  with  cell-to-cell  variation 
can  be  represented  as  a  linear  combination  of  the  mode  shapes  of  the  nominal  system,  and  (3)  the  frame 
holding  each  cell  has  vibratory  motion.  A  numerical  example  of  an  academic  battery  pack  with  pouch 
cells  is  presented  to  demonstrate  that  the  PROM  captures  the  effects  of  both  prestress  and  structural 
variation  on  battery  packs.  The  PROM  is  validated  numerically  by  comparing  full-order  finite  element 
models  (FEMs)  of  the  same  systems. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Typical  hybrid  electric  vehicle  (HEV)  battery  packs  are  assem¬ 
bled  by  bolts  or  welds  for  joining  cells  within  the  pack  and  for 
integrating  the  battery  structure  into  the  rest  of  the  vehicle.  The 
prestress  due  to  joining  can  affect  the  dynamic  response  of  the 
structures  significantly.  In  addition,  HEV  battery  pack  structures 
typically  include  dozens  or  even  hundreds  of  cells.  Because  the  cells 
are  nominally  identical,  if  the  cells  are  arranged  in  a  spatially 
repeating  layout  then  the  battery  pack  falls  under  the  class  of 
structures  known  as  periodic  structures.  The  dynamics  of  periodic 
structures  are  known  to  feature  very  high  modal  density  in  many 
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frequency  bands.  The  high  modal  density  combined  with  small, 
random  structural  variations  among  the  cells  (which  are  unavoid¬ 
able  in  practice)  can  lead  to  drastic  consequences  on  the  structural 
dynamics.  In  particular,  they  can  exhibit  Anderson  localization  [  1  — 
3  ],  where  the  vibration  energy  is  spatially  confined  to  a  small  region 
of  the  structure. 

The  focus  of  this  work  is  to  develop  computationally  efficient 
methods  for  predicting  the  vibratory  response  of  a  battery  pack  that 
may  suffer  from  localization.  It  is  believed  that  certain  battery 
packs,  such  as  those  with  pouch  cells,  may  be  especially  susceptible 
to  localization  phenomena.  Localization  can  lead  to  battery  cell 
damage  due  to  sharp  increases  in  local  vibration  amplitude  and 
stress  levels  in  a  pack,  which  in  turn  can  lead  to  mechanical  or 
electromechanical  failure. 

Stress  and  failure  issues  in  battery  packs  have  been  investigated 
in  several  previous  studies.  For  example,  Offer  et  al.  [4]  provided  a 
new  concept  of  module  design  and  a  fault  diagnosis  technique  in 
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electric  vehicle  batteries.  They  focused  on  testing  the  full  battery 
pack  and  diagnosing  subsequent  problems  related  to  cells  to 
demonstrate  how  a  full  vehicle  test  can  be  used  to  identify  mal¬ 
functioning  strings  of  cells.  Xiao  et  al.  [5]  developed  a  multi-scale 
approach  for  the  stress  analysis  of  polymeric  separators  in  a 
lithium-ion  battery.  They  presented  a  finite  element  method  based 
on  a  multi-scale  approach  for  the  stress  analysis  of  the  separator  in 
a  battery  cell.  Sahraei  et  al.  [6]  revealed  the  underlying  physics  and 
identified  important  groups  of  parameters.  To  do  that,  they  per¬ 
formed  tests  on  pouched  and  bare  lithium-ion  cells  under  several 
loading  conditions.  They  validated  the  test  by  numerical  simula¬ 
tions.  The  outcomes  from  the  test  data  were  used  for  the  devel¬ 
opment  of  an  advanced  constitutive  model  needed  for  strength/ 
weight  optimization  and  safety  assessments  of  Li-ion  batteries. 
Kenney  et  al.  [7]  established  the  relationship  between  cell-to-cell 
variations  due  to  their  manufacturing  process  and  their  overall 
impact  on  a  full  battery  pack.  They  used  mathematical  modeling 
techniques  to  predict  how  slight  variations  in  the  manufacturing  of 
electrodes  manifest  themselves  in  an  operational  battery  module. 
The  variation  in  the  performance  of  the  module  has  been  quantified 
as  a  function  of  the  manufacturing  variation  at  the  electrode  level. 
Dubarry  et  al.  [8  j  developed  modeling  and  simulation  techniques  to 
show  that  accurate  battery  pack  simulations  can  be  achieved  if  cell- 
to-cell  variations  are  taken  into  account.  They  used  an  accurate 
single-cell  model  which  was  validated  against  experimental  data  to 
provide  high  fidelity  pack  simulation. 

However,  the  structural  characteristics  of  the  entire  battery  pack 
with  very  high  modal  density  have  not  been  previously  explored. 
Therefore,  the  literature  does  not  address  the  possible  occurrence 
or  prediction  of  localization  in  battery  packs  due  to  random 
structural  variations  among  the  cells. 

In  order  to  properly  account  for  the  effects  of  random  variations, 
statistical  dynamic  response  calculations  are  needed.  Such  statis¬ 
tical  calculations  are  hard  to  perform  using  linear  methods  because 
the  mode  shapes  of  a  pack  depend  in  a  nonlinear  fashion  on  the 
parameters  of  each  cell.  The  alternative  is  to  use  sample-based 
statistical  analyses,  such  as  Monte  Carlo  simulations.  However, 
typical  finite  element  models  (FEMs)  of  battery  packs  have  several 
million  degrees  of  freedom  (DOF).  Thus,  the  computational  time  for 
obtaining  just  a  single  sample  can  be  on  the  order  of  a  day.  This  may 
make  it  infeasible  to  perform  Monte  Carlo  simulations  with  FEMs. 

To  overcome  this  issue,  in  the  field  of  structural  dynamic  anal¬ 
ysis,  component  mode  synthesis  (CMS)  [9-15]  is  well  established 
as  an  alternative  to  conventional  FEMs  with  large  numbers  of  DOF. 
CMS  belongs  to  a  wide  class  of  domain  decomposition  techniques. 
In  a  CMS  approach,  the  global  structure  is  divided  into  several 
substructures,  and  the  DOF  of  each  individual  substructure  are 
reduced  significantly.  Then,  the  substructures  are  reconnected,  and 
the  dynamic  response  of  the  system  is  predicted  very  efficiently 
and  accurately.  However,  classical  CMS  must  be  modified  in  order 
to  account  for  parametric  variability  in  the  structure.  Thus,  alter¬ 
nate,  design-oriented  techniques  have  been  developed.  One  such 
approach  is  to  generate  what  is  referred  to  as  parametric  reduced- 
order  models  (PROMs).  PROMs  were  introduced  initially  by  Balmes 
et  al.  [16,17]  to  avoid  the  expensive  process  of  reanalysis  of  complex 
structures.  In  addition,  several  other  PROM  methods  have  been 
developed  [18-22], 

In  particular,  the  multiple-component  PROMs  (MC-PROMs)  [21  ] 
have  been  developed  by  the  authors.  MC-PROMs  are  well  suited  for 
the  structure  modeled  with  shell-type  finite  elements.  However, 
MC-PROMs  have  several  drawbacks,  namely:  (1)  a  numerical 
instability  of  system  level  matrices  can  be  encountered  due  to  the 
transformation  matrix,  (2)  MC-PROMs  do  not  capture  well 
elemental  level  nonlinearities  for  thickness  variations  of  brick  type 
finite  elements,  and  (3)  the  interface  DOF  are  hard  to  reduce.  Thus, 


the  next-generation  PROMs  (NX-PROMs)  technique  have  been 
developed  recently  to  overcome  these  drawbacks  [22],  The  con¬ 
cepts  used  in  NX-PROMs  are  applied  herein  to  capture  the  prestress 
variations. 

The  PROM  techniques  are  highly  efficient  methods  for  esti¬ 
mating  the  statistics  of  the  structural  dynamic  response.  However, 
for  a  structure  with  very  high  modal  density,  previously  developed 
PROMs  have  to  be  modified  to  efficiently  capture  the  dynamic 
response.  In  particular,  such  a  modification  can  be  similar  to  the 
component  mode  mistuning  (CMM)  [23]  method.  CMM  was 
developed  for  predicting  the  dynamic  response  of  bladed  disks 
found  in  turbo-machinery.  Typically,  these  bladed  disks  suffer  from 
high  modal  density.  Thus,  small  structural  variations  in  the  blades 
significantly  affect  the  system-level  dynamic  response.  Nonethe¬ 
less,  it  has  been  shown  that  the  mode  shapes  of  a  mistuned  bladed 
disk  can  be  represented  as  a  linear  combination  of  the  mode  shapes 
of  the  tuned  bladed  disk  [24],  This  allows  CMM  to  capture  the 
dynamic  response  effectively  with  a  small  number  of  DOF,  and  that 
is  the  inspiration  for  the  new  PROMs  introduced  in  this  work. 

This  paper  is  organized  as  follows.  In  Section  2,  structural 
properties  and  dynamic  characteristics  of  HEV  battery  packs  are 
evaluated.  Next,  in  Section  3,  PROM  is  applied  for  prestress  vari¬ 
ations,  and  new  concept  of  PROM  is  developed  for  cell-to-cell 
structural  variations.  In  Section  4,  numerical  examples  of  an  aca¬ 
demic  battery  pack  models  are  used  to  demonstrate  and  test  the 
proposed  methods,  including  a  Monte  Carlo  simulation  technique 
that  supports  fatigue  life  prediction  by  quantifying  stochastic  ef¬ 
fects  in  the  forced  response  and  sensitivities  to  various  sources  of 
parameter  variation.  Finally,  conclusions  are  summarized  in  Sec- 


2.  Structural  properties  and  dynamic  characteristics  of 
battery  packs 

HEV  battery  packs  typically  have  100—300  individual  cells  that 
are  nominally  identical.  To  demonstrate  the  structural  dynamic 
characteristics  of  a  battery  pack,  a  simplified  academic  model  of  a 
pack  of  pouch  cells  was  developed  using  finite  elements  as  shown 
in  Fig.  1.  The  nodes  on  one  end  of  the  pack  are  totally  fixed.  The 
prestress  and  dynamic  loads  are  applied  to  nodes  on  the  other  end 
of  the  structure.  The  prestress  loads  are  applied  to  the  longitudinal 
direction  to  compress  the  structure,  and  the  harmonic  excitations 
are  applied  to  all  three  (x,  y,  z)  directions.  The  excitation  frequencies 
are  in  the  range  of  2200—2350  Hz.  A  total  of  20  nominally  identical 
cells  are  stacked.  Typically,  foam  or  epoxy  stiffness  layers  can  be 
placed  between  cells.  These  elements  are  nonlinear  materials,  but 
they  are  soft.  Thus,  herein  these  nonlinear  materials  are  ignored. 

To  evaluate  how  the  prestress  variations  in  the  structure  affect 
the  structural  response,  two  different  levels  of  prestress  were 
applied:  150  kN  and  470  kN.  Fig.  2  shows  the  forced  responses  of 
the  center  node  of  the  1st  cell  for  the  different  levels  of  prestress. 
The  responses  obtained  by  full-order  models  are  significantly 
different  for  each  level  of  prestress.  For  the  prestress  levels,  the 
intent  is  to  demonstrate  that  the  PROMs  can  capture  the  structural 
dynamics  for  different  prestress  values  that  lead  to  a  significant 
change  in  the  forced  response.  The  values  150  kN  and  470  kN  were 
chosen  because  they  met  these  criteria  for  demonstrating  the 
PROM  performance.  However,  it  should  be  noted  that  these  are  just 
two  sampled  prestress  values,  and  the  precise  numbers  do  not  have 
special  significance. 

Fig.  3  shows  a  single  cell  and  the  frames  that  join  it  to  the 
adjacent  cells.  The  single  pouch  cell  is  a  plate-like  structure  for 
which  the  small,  random  structural  variations  are  modeled  by 
perturbations  to  the  elastic  modulus  (£).  The  (nominally  identical) 
cells  are  mechanically  coupled  through  the  frames,  which  induces  a 
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Fig.  1.  The  geometry  of  a  battery  pack  with  20  pouch  ce 
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high  modal  density  to  the  entire  battery  pack  structure,  as  shown  in 
Fig.  4.  The  flat  regions  in  Fig.  4  indicate  frequency  ranges  of  high 
modal  density.  For  example,  there  are  over  20  modes  in  the  range 
2200-2300  Hz.  If  the  battery  pack  had  more  cells,  the  modal 
density  would  be  even  higher.  The  frequency  range  2200-2300  Hz 
features  high  modal  density,  as  shown  in  Fig.  4.  In  such  a  frequency 
range,  the  forced  response  of  a  nominally  periodic  structure  can  be 
sensitive  to  small,  random  variations  [2,3],  That  is  the  reason  for 
setting  the  excitation  frequency  in  the  range  of  2200-2300  Hz. 
Note  that  there  is  also  flat  region  of  natural  frequencies  about 
4000  Hz,  and  one  would  also  expect  to  observe  sensitivity  of  the 
forced  response  to  parameter  variation  in  this  frequency  range. 
However,  for  this  academic  example,  there  was  no  specific  reason 
to  examine  a  higher  frequency  range.  Therefore,  the  lower  fre¬ 
quency  range  was  selected  for  this  initial  study. 

To  examine  how  structural  variations  in  the  cells  affect  the 
structural  response,  Young’s  modulus  variations  were  applied,  as 


Frequency  (Hz) 


Upper 


Single  battery  cell 
(structural  variation) 


Fig.  3.  The  geometry  of  a  single  pouch  cell. 


described  in  Table  1.  Young’s  modulus  changes  were  used  as  a 
simple  way  to  vary  the  cell  stiffnesses  directly  in  both  the  finite 
element  model  and  the  PROM.  In  practice,  some  small  differences 
among  the  stiffnesses  of  the  cells  are  unavoidable  due  to  thickness 
variations  and  other  manufacturing  tolerances.  In  addition,  even  if 
each  cell  started  perfectly  identical,  as  the  battery  cells  are  charged 
and  discharged  they  experience  large  temperature  fluctuations, 
and  the  individual  cell  stiffnesses  will  vary  accordingly  during 
operation.  Therefore,  Young’s  modulus  changes  were  used  to  model 
the  overall  variability  in  cell  stiffnesses.  Furthermore,  in  order  to 
observe  the  consequences  of  small  parameter  variations  in  just  a 
few  cells,  for  cases  1  and  2  the  variations  were  applied  in  only  four 
cells.  However,  for  other  cases,  the  variations  were  applied  to  all 
cells,  which  is  the  more  realistic  case. 

The  mode  shapes  of  the  structure  with  nominal  parameters  (no 
variation)  were  compared  with  the  mode  shapes  of  the  structure 
for  cases  1  and  2.  These  mode  shapes  are  shown  in  Fig.  5.  Although 
the  structural  variations  were  small,  the  mode  shapes  are  affected 
significantly.  In  particular,  note  that  some  mode  shapes  are  local¬ 
ized  about  a  few  cells  in  which  there  are  no  parameter  variations. 

The  forced  response  for  the  nominal  case  was  also  compared  to 
the  forced  response  for  each  case  of  variation,  as  shown  in  Fig.  6. 
These  results  show  that  small  local  structural  variations  can  induce 
large  changes  in  the  global  forced  response.  Fig.  6  shows  the 
response  of  the  4th  cell  (top  plot)  and  the  6th  cell  (bottom  plot).  As 
shown  in  Table  1,  there  are  no  variations  in  the  parameters  of  cells  4 


Fig.  2.  Forced  response  of  the  center  node  of  the  1  st  cell  for  different  levels  of  prestress  Fig.  4.  Natural  frequencies  of  the  academic  battery  pack  without  prestress  and  without 
without  structural  variation.  cell-to-cell  variations. 
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Table  1 

Two  patterns  of  Young's  modulus  variations. 

Pattern  1  Pattern  2 

Cell  Variation  Cell  Variation 

15%  33% 

5  -7%  9  -5% 

12  1%  13  2% 

16  3%  20  -5% 


and  6.  However,  there  are  significant  changes  in  the  dynamic 
response  of  these  cells.  For  example,  the  maximum  response  of  the 
6th  cell  is  approximately  0.0085  mm,  whereas  the  maximum 
response  of  the  same  cell  in  case  2  of  parameter  variations  is 
approximately  0.0061  mm.  This  is  almost  a  30%  increase  in  the 
maximum  response,  even  though  the  maximum  parameter  varia¬ 
tion  between  the  nominal  battery  and  that  of  case  2  is  only  5%.  This 
demonstrates  that  small  local  structural  variations  can  have  large 
global  consequences.  To  efficiently  capture  the  dynamic  charac¬ 
teristics  of  the  structure  with  simultaneous  prestress  and  structural 
variations,  a  new  approach  was  developed,  as  described  in  the  next 
section. 

3.  Parametric  reduced-order  model 

3.1.  Modeling  prestress  in  a  structure 

The  prestress  affects  the  dynamic  response  of  structures  by 
changing  their  stiffness.  The  first  key  idea  of  PROMs  is  the  param¬ 
eterization  of  the  stiffness  matrices  [22],  The  parameterization 
significantly  reduces  the  reanalysis  time  because  any  variation  in 


the  parametric  range  can  be  applied  for  different  levels  of  prestress 
without  full-order  finite  element  calculations  or  mesh  refinement. 
The  other  key  idea  is  to  predict  vibration  responses  using  reduced- 
order  models  (ROMs)  as  opposed  to  full-order  models  to  reduce  the 
calculation  time.  To  detail  the  construction  of  PROMs,  these  two  key 
ideas  are  implemented. 

The  transformation  matrix  for  PROMs  is  constructed  by  a  set  of 
fixed-interface  normal  modes  <l>aUg  given  by 

<S  =  [*0  <  *2  <:> 

where  matrix  <l>n  corresponds  to  the  nominal  parameter  values, 
and  matrices  (i  =  1,  2,  3)  correspond  to  3  other  parameter 
values.  A  third-order  interpolation  is  used  for  a  parameterization. 
Thus,  4  samples  of  matrices  (at  distinct  parameter  values)  are 
needed  for  constructing  the  transformation  matrix. 

In  general,  the  columns  of  4>jUg  are  not  orthogonal.  Therefore, 
for  numerical  stability,  an  orthogonal  basis  for  the  space  spanned 
by  these  modes  is  computed.  This  basis  is  obtained  by  a  truncated 
set  of  left  singular  vector  UN  of  <l>aUg  [21,22],  The  basis  of  vectors  in 
UN  is  orthogonal,  so  numerical  conditioning  problems  are  avoided. 
In  addition,  the  left  singular  vectors  UN  are  truncated  to  lower  the 
size  of  the  resulting  model.  The  choice  of  the  cutoff  point  can  affect 
accuracy  if  it  is  chosen  too  high,  but  it  does  not  affect  numerical 
stability.  Thus,  the  cutoff  value  can  be  estimated  by  a  standard 
convergence  study  where  calculations  are  done  at  ever  lower  cutoff 
values  until  convergence  is  obtained.  For  the  calculations  per¬ 
formed  in  this  study,  the  cutoff  was  0.01%  of  the  maximum  singular 
value.  Then  UN  is  used  in  a  transformation  matrix  tprom  to  convert 
for  physical  coordinates  to  generalized  coordinates.  The  trans¬ 
formed  stiffness  matrices  for  the  4  samples  of  such  matrices 


(a)  Nominal 


(b)  Case  1 


(c)  Case  2 


Fig.  5.  Mode  shapes  of  the  nominal  structure  with  identical  cells  (left)  and  the  structures  with  cases  1  (center)  and  2  (right)  of  cell  parameter  variations. 
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Fig.  6.  Forced  response  of  the  center  node  of  the  4th  cell  (top)  and  the  6th  cell  (bot¬ 
tom)  for  the  nominal  case  of  no  parameter  variations  versus  two  cases  with  parameter 
variations. 


rKo.eql 

'  K(p0)eq  ‘ 

Kl.eq 

K(p0  +  5p)eq 

K2,eq 

I<(P0  +  2<5p)eq 

_  K3,  eq 

K(p0+3<5p)eq 

where 


10  0  0 

1  5p  (dp)2  (dp)3 

1  25 p  (25 p)2  (25 p)3 
.1  35 p  (35p)2  (35 p)3. 


Equation  (4)  can  be  easily  solved  by  simply  inverting  the  4x4 
matrix  C.  Let  us  denote  by  A  this  inverse  matrix,  i.e. 


A  =  C  1 


An  A12  A13  A14 
A21  A22  A23  A24 
A3i  A32  A33  A34 

_a41  a42  a43  a44 


Rearranging  Eq.  (4)  using  the  entries  A, 


obtains 


K(p0  +  Ap)~K™OM 

=  b0lC°M  +  b]  K|JR0M  +  fa2K2ROM  +  b3KPR0M, 

where 

bo  =  (An  +  A21Ap+A31Ap2+A41Ap3), 
b\  =  (A12+A22Ap+A32Ap2+A42Ap3), 
b2  =  (An  +  A23  A p  +  A33Ap2  +  A43Ap3) , 

=  (Ai4  +  A24Ap  +  A34Ap2  +  A44Ap3) . 

The  parameterized  stiffness  matrix  in  Eq.  (5)  is  used  to  capture 
the  dynamic  response  of  the  structure  with  different  levels  of 
prestress.  Details  of  this  procedure  can  be  found  also  in  Hong  et  al. 
[22], 

This  section  described  the  way  to  capture  the  prestress  effect 
efficiently  based  on  NX-PROMs  [22],  In  the  next  section,  the  nu¬ 
merical  methods  to  predict  the  dynamic  response  of  the  periodic 
structure  which  features  very  high  modal  density  are  formulated. 


kPROM  _  ^tPROM^TkFEMjPROM 

where  (i  =  0, 1,  2, 3). 

For  convenience,  the  parameterization  procedure  is  summa¬ 
rized  next.  This  procedure  closely  follows  that  of  Hong  et  al.  [22], 
The  parameterization  equation  consists  of  a  third-order  interpola¬ 
tion,  which  can  be  written  as 

K(p0  +  Ap)  «K0  +  ApK,  +  Ap2K2  +  Ap3I<3. 

Four  equations  are  needed  to  calculate  the  matrices  K,  (for  i  =  0, 
1,  2,  3).  To  obtain  these  four  equations,  reduced-order  stiffness 
matrices  are  computed  using  Eq.  (1)  for  four  parameter  values. 
First,  consider  the  case  where  A p  =  0.  One  obtains 

K(p0)  ~KoROM-  (2) 

Next,  consider  Ap  =  idp  (with  i  =  1,  2,  3),  one  obtains 

K(Po  +  i5p) » I<rR0M  =  K0  +  (tfp)Kt  +  (i<5p)2K2  +  (i5p)3K3 .  (3 ) 

Rearranging  Eqs.  (2)  and  (3)  into  matrix  form,  for  each  entry  e,  q 
of  the  matrices  Ko.  Ki,  K2  and  K3,  one  obtains 


3.2.  Modeling  structural  variations  in  cells 

The  equations  of  motion  for  the  structure  with  no  variation  and 
with  variation  can  be  expressed  as 

Mx  +  Cx  +  Kx  =  F,  (6) 

Mx+Cx+  (K+K^x  =  F,  (7) 

where  K'5  contains  the  stiffness  variation.  Based  on  Eqs.  (6)  and  (7), 
the  mode  shapes  are  defined  by  the  following  eigenvalue  problems 

K &  -  =  0, 

(K  +  K*)$m  -  =  0, 

where  superscript  t  and  m  indicate  the  tuned  (nominal)  and  mis- 
tuned  (structure  with  variation)  quantities.  Matrices  Ar  and  Am  are 
diagonal  and  contain  the  eigenvalues  of  the  tuned  and  mistuned 
systems.  If  the  structure  has  mass  variations,  the  mass  matrix  in  Eq. 
(7)  is  (M  +  M5).  For  the  sake  of  simplicity,  only  the  case  of  stiffness 
variations  was  considered. 

The  novel  approach  is  based  on  two  key  assumptions.  The  first 
assumption  is  that  the  mode  shapes  4>m  of  a  pack  with  parametric 
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Fig.  7.  1st  and  2nd  modes  of  a  fixed-boundary  pouch  cell. 


variations  can  be  approximated  as  a  linear  combination  of  the 
mode  shapes  <I>t  of  nominal  pack  with  no  parametric  variations. 
This  first  assumption  is  ensured  by  the  high  modal  density. 

The  second  assumption  for  structural  variations  is  that  the 
variations  in  stiffness  of  a  cell  can  be  projected  onto  a  small  set  of 
modes  of  the  nominal  cell  with  a  fixed  boundary,  as  shown  in  Fig.  7. 
This  second  assumption  relies  on  the  fact  that  the  boundary  mo¬ 
tions  can  be  ignored  when  the  boundary  of  the  pouch  cell  is  not 
moving.  However,  it  turns  out  that  the  boundary  motion  has  to  be 
considered  because  it  is  not  small,  as  shown  in  Fig.  8.  Thus,  the 
plate-like  modes  of  a  nominal  cell-with  its  boundary  displaced  the 
same  amount  as  the  frame-are  used  in  the  proposed  PROMs.  This  is 
a  key  step  for  ensuring  accuracy.  This  approach  is  distinct  from  the 
original  CMM  method  [23]  because  the  CMM  method  does  not 
account  for  the  boundary  motion. 

These  two  key  ideas  are  implemented  to  model  a  system  with 
prestress  variations.  Eq.  (6)  for  the  tuned  (nominal)  structure  with 
prestress  variations  can  be  written  as 

M™°0muprom  +  (1  +  jT)K^°oMuPROM  =  Fprom,  (8) 

where  uPR0M  are  generalized  coordinates  given  by  x  =  tpromuprom. 
Subscript  A p  indicates  the  prestress  variation  and  subscript  0  in¬ 
dicates  no  structural  variations.  7  is  a  structural  damping 
coefficient. 

The  eigenvalue  problem  determined  from  Eq.  (8)  with  fprom  =  0 
is  solved  to  determine  the  natural  frequencies  APR°M  and  mode 


shapes  <f>^p°M  of  the  system.  Next,  a  secondary  modal  analysis  can 
be  used  to  transform  coordinates  from  uPR0M  to  vPR0M  such  that 
uprom  _  *promvpr0m  Eq  (8)  becomes 

vPROM  +  (l+j7)APRO0MvPROM  =  ($PRO0M)VROM. 

Note  that  a£r°m  contains  eigenvalues  of  the  prestressed  struc¬ 
ture  without  variation.  That  is  different  from  Af  which  contains 
eigenvalues  of  the  tuned  system  without  prestress. 

As  described  above,  the  first  key  assumption  is  used  in  the  PROM 
domain.  The  mistuned  mode  shapes  of  the  structure  can  be  rep¬ 
resented  as  a  linear  combination  of  tuned  mode  shapes.  Thus,  one 
can  assume  the  physical  coordinates  x  in  Eq.  (7)  can  be  expressed  as 

x  =  jPROMyPROM  yPROM<j,PROMvPROM 

Then  the  equations  of  motion  of  the  mistuned  structure  in  the 
vprom  coorcjinates  can  be  written  as 


vPROM  +  (l+j7)(APR00MvPROM+A)  =  (^°0M)TFprom, 

where  A  =  (^pr7)t(Tprom)tKpe^Tprom4»pro0mvprom,  and  Kp™ 
contains  the  variations  in  the  stiffness  matrix  in  physical  co¬ 
ordinates  between  the  nominal  structure  and  the  structure  with 
variations. 

To  parameterize  the  stiffness  variation  in  the  jth  cell,  a  first- 
order  interpolation  was  used,  which  is  given  by 

=  "/O 

where  J  is  the  jth  cell  stiffness  matrix  for  the  nominal  cell,  and 
rrij  is  the  amount  of  stiffness  variations  in  the  jth  cell.  Then  K™J 
can  be  partitioned  into 


\1<°M  <0 

[  koa  koo  J 


where  subscript  A  and  0  are  DOF  refer  to  interface  (A)  and  internal 
(0)  DOF  of  the  jth  cell.  The  partitions  of  the  stiffness  matrix  are  used 
to  calculate  the  static  deformation  from  the  boundary  displacement 
of  the  jth  cell  by  using 


where  T,,  and  T0,  are  interface  and  internal  displacements  of  the 
jth  cell.  The  interface  displacement  T/y  is  obtained  from  the 
transformation  matrix  TPR0M  by  selecting  entries  corresponding  to 
interface  DOF  of  the  jth  cell  from  global  DOF,  which  is 


Taj  =  T™0M . 

Then  the  static  deformation  induced  by  the  boundary 
displacement  can  be  computed  by  Eq.  (9),  that  is 


[jPROM  I 

-K°00’CTPrJ 


Next,  a  set  of  fixed  interface  normal  modes  <FC  of  a  single  battery 
cell  structure  are  calculated.  The  fixed-interface  normal  modes  are 
a  truncated  set  of  modes  obtained  by  solving  the  eigenproblem 
with  the  mass  and  stiffness  matrices  of  the  cell  with  a  fixed 
boundary.  One  obtains 
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Koodoo  ~~  <0*00^  —  0) 
where  <f>c  =  [$m/$So]  =  [°/*oo]  ■ 

Next,  the  computed  static  deformation  T,  and  a  set  of  fixed 
interface  normal  modes  4>c  of  the  single  cell  are  used  to  calculate 
modal  participation  factors  q,  as  follows 


These  relations  are  used  to  construct  the  equations  of  motion  for 
prestress  and  stiffness  variations  as  follows 

vPR0M  +  (1  +  Jy)  Apro0mvprom  +  B  =  (4>^°0m)TFprom, 

where 

B  =  (1  +  jy)  (^p°oM)TC^^°oMvPROM, 

c=  \jr  ($cqj  +T,-)TKPEp7'(«hcqj + Tj)| , 

and  N  is  the  number  of  cells  that  have  stiffness  variations. 

The  PROM  for  the  HEV  battery  pack  structure  has  been  formu¬ 
lated  to  predict  the  dynamic  response  when  the  structure  has  a 
prestress  and  cell-to-cell  variations.  In  the  next  section,  the 
computational  efficiency  and  accuracy  of  the  newly  developed 
PROM  is  evaluated  by  applying  it  to  the  academic  battery  pack. 


4.  Numerical  results 

As  a  benchmark  for  testing  the  performance  of  the  proposed 
method,  forced  response  results  were  been  obtained  using  the 
academic  battery  finite  element  model  shown  in  Fig.  1.  This  FEM 
has  208,753  DOF  and  features  20  nominally  identical  cells. 

The  frequency  range  of  interest  is  1500-3000  Hz,  which  is  the 
first  flat  region  of  the  natural  frequencies  plotted  in  Fig.  4.  This  type 
of  frequency  band  with  high  modal  density  is  characteristic  of  pe¬ 
riodic  structures.  For  system  modes  in  such  a  region,  the  vibration 
response  shape  of  each  cell  will  tend  to  be  very  close  to  that  of  a 
mode  shape  for  a  single,  isolated  cell.  For  this  particular  band  of 
system  modes,  the  response  shape  of  each  cell  is  very  close  to  the 
first  fixed-boundary  mode  shape  of  a  cell,  which  is  the  top  mode 
shape  shown  in  Fig.  7. 

Prestress  and  dynamic  loads  were  applied  as  shown  in  Fig.  1.  For 
the  dynamic  loads,  harmonic  forcing  was  applied  at  the  right  end 
(cell  20)  while  the  boundary  DOF  at  the  left  end  (cell  1)  were  held 
fixed.  For  this  initial  study,  both  the  model  and  the  loading  were 
intentionally  kept  simple.  More  realistic  models  and  loading  con¬ 
ditions  will  be  examined  in  future  work. 

For  the  following  numerical  results,  the  forced  response  ampli¬ 
tudes  were  obtained  at  the  center  nodes  of  the  cells.  As  with  the  free 
response,  the  forced  response  shape  of  each  cell  was  very  close  to  the 
first  mode  shape  of  an  individual  cell  (multiplied  by  a  scalar 
amplitude),  based  on  direct  observation  of  the  finite  element  anal¬ 
ysis  results.  In  particular,  the  response  shape  of  each  cell  was  similar 
to  that  of  the  first  fixed-boundary  mode  shape.  (It  is  recognized  that 
this  might  lead  one  to  assume  that  each  cell  will  not  interact 
dynamically  with  neighboring  cells.  However,  it  should  be  noted 


Table  2 

Comparison  of  the  full-order  models  and  the  PROM. 

Full-order  model  PROM 

System  DOF  208,753  38 

Reanalysis  time  10,120-10,780  s  0.89-0.95  s 


that  the  boundary — the  frame — is  not  fixed.  In  fact,  the  frame  vi¬ 
brates  and  provides  significant  dynamic  coupling  among  the  cells.) 
Therefore,  the  maximum  forced  response  displacement  for  each  cell 
occurred  at  the  center  node.  Furthermore,  the  amplitudes  at  the 
center  nodes  of  the  cells  are  representative  of  not  only  the  maximum 
displacement  levels  but  also  the  maximum  stress  levels  of  the  cells. 
Therefore,  the  effects  of  prestress  and  structural  variation  on  these 
amplitudes  provide  an  excellent  indicator  of  the  effects  on  fatigue 
life.  Higher  amplitudes  indicate  correspondingly  higher  stresses  as 
well  as  an  increased  risk  of  structural  fatigue  failure. 


4.1  Validation  of  PROM  performance 

For  simulating  the  forced  response  in  the  frequency  range  of 
interest,  a  38-DOF  PROM  was  generated  from  the  full-order  FEM 
(which  has  208,753  DOF).  The  number  of  DOF  and  the 


Case  1  for  150kN 


cell  (bottom)  predicted  by  the  FEM  and  a  PROM  for  cases  1  (top)  and  2  (bottom)  with 
150  kN  prestress. 
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computational  time  required  for  the  FEM  and  PROM  are  shown  in 
Table  2.  The  reanalysis  time  for  each  variation  is  about  10,000  times 
shorter  for  the  PROM  compared  to  the  FEM. 

This  computational  gain  would  increase  if  a  finer  finite  element 
mesh  were  used.  That  is  because  the  PROM  captures  the  low¬ 
dimensional  physics  of  the  problem.  This  low  dimensionality 
means  that  relatively  few  coordinates  are  necessary  to  describe  the 
dynamics  of  the  actual  physical  system.  This  number  of  coordinates 
is  a  feature  of  the  physics,  not  of  the  model  used  to  discretize  the 
physics.  The  finite  element  model  will  increase  in  size  as  the  mesh 
is  refined;  however,  the  physics  remains  the  same  and  are  captured 
by  the  same  PROM  coordinates.  Therefore,  the  computational  sav¬ 
ings  from  using  the  PROM  are  expected  to  be  much  greater  for  a 
more  realistic  FEM  of  an  industrial  battery  pack. 

As  example  cases  for  validating  the  accuracy  of  the  PROM,  two 
patterns  of  Young’s  modulus  variations  were  applied  to  the  cells,  as 
shown  in  Table  1.  In  addition,  for  each  variation  pattern,  prestress 
values  of  150  kN  and  470  kN  were  applied.  Forced  response  simu¬ 
lations  were  then  run  by  sweeping  the  excitation  through  the  fre¬ 
quency  range  of  interest  for  all  four  combinations  of  prestress 
values  and  variation  patterns. 

Figs.  9  and  10  shows  the  response  of  the  center  node  of  the  5  th  cell 
for  variation  pattern  1,  and  the  9th  cell  for  variation  pattern  2  under 
both  150  kN  and  470  kN  prestress.  The  solid  lines  indicate  predictions 


Case  1  for  470kN 


9th  cell  (bottom)  predicted  by  the  FEM  and  a  PROM  for  cases  1  (top)  and  2  (bottom) 
with  470  kN  prestress. 


of  full-order  models  and  the  symbols  indicate  PROM  predictions.  For 
all  of  these  forced  response  simulations,  it  can  be  seen  that  the  PROM 
results  agreed  very  well  with  the  much  more  expensive  FEM  results. 
Figs.  11  and  12  show  the  maximum  error  between  the  PROM  and  the 
FEM  predictions  for  all  20  cells  over  the  entire  frequency  range  of 
interest.  The  maximum  errors  were  no  larger  than  2.7%,  and  most 
were  below  1%  for  all  four  simulation  cases. 

4.2.  Effects  of  prestress  for  nominal  system 

Fig.  13  shows  the  maximum  forced  response  amplitude  in  the 
frequency  range  1500-3000  Hz  for  each  cell  in  the  nominal  system. 
It  is  important  to  note  that  these  amplitudes  are  the  maximum 
values  for  the  1500—3000  Hz  range,  and  the  cells  do  not  necessarily 
experience  these  amplitudes  at  the  same  excitation  frequencies. 
The  contributions  from  multiple  system  resonances  may  be  rep¬ 
resented  in  this  plot.  That  is,  these  plots  are  not  resonant  response 
shapes.  Rather,  they  are  the  baseline  values  for  the  maximum  vi¬ 
bration  response  experienced  by  each  cell  in  the  nominal  battery 
pack.  They  serve  as  a  baseline  for  in  the  absence  of  cell-to-cell 
structural  variations. 
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Cell  number 
Case  2  for  470kN 
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Cell  number 


Fig.  12.  Maximum  errors  between  PROM  and  FEM  for  each  cell  (center  node 
displacement)  in  the  frequency  range  of  interest  for  cases  1  (top)  and  2  (bottom)  with 
470  kN  prestress. 


These  forced  response  simulation  results  are  shown  in  Fig.  13  for 
two  different  prestress  levels,  150  kN  and  470  kN.  It  can  be  seen  that 
the  maximum  amplitudes  versus  cell  number  show  a  similar  trend. 
As  expected,  for  both  cases,  the  highest-responding  cells  are  those 
closest  to  the  excitation  source.  However,  the  different  prestress 
levels  lead  to  noticeable  differences  in  the  maximum  amplitudes 
for  these  highest-responding  cells. 

It  should  be  noted  that  the  prestress  level  influences  the 
strength  of  the  structural  coupling  between  cells.  It  is  well  known 
that  the  amount  of  Anderson  localization  in  nearly  periodic  struc¬ 
tures  depends  not  only  on  the  amount  of  structural  variations  but 
also  on  the  strength  of  coupling  [2,3],  Therefore,  one  would  also 
expect  to  see  different  effects  of  structural  variations  for  different 
prestress  levels.  This  will  be  examined  next. 


4.3.  Effects  of  prestress  and  structural  variations 

In  a  battery  pack,  a  single  cell  with  intense  vibrations  may  lead 
to  the  failure  of  the  entire  pack.  Thus,  the  maximum  response  level 
for  any  cell  in  a  battery  pack  is  a  key  reliability  metric.  However,  in 


Forced  Response  of  Nominal  System  (Identical  Cells) 


Fig.  13.  Maximum  forced  response  amplitudes  for  the  nominal  system. 


order  to  consider  the  effects  of  random  structural  variations,  a 
probabilistic  approach  is  needed.  For  example,  a  Monte  Carlo 
simulation  with  thousands  of  variation  patterns  might  be  required. 
For  a  large  FEM  of  an  industrial  battery  pack,  Monte  Carlo  simula¬ 
tions  may  be  prohibitively  expensive.  In  contrast,  the  accuracy  and 
high  computational  speed  of  the  PROM  make  it  feasible  to  run 
Monte  Carlo  simulations. 

In  this  study,  Monte  Carlo  simulations  were  performed  using  the 
38-DOF  PROM  to  calculate  the  forced  response  for  10,000  different 
patterns  of  cell-to-cell  elastic  modulus  variations.  The  variation 
patterns  were  sampled  from  a  normal  distribution  with  zero  mean 
and  standard  deviation  set  to  5%  of  the  nominal  value.  This  corre¬ 
sponds  to  a  standard  deviation  of  approximately  2.5%  with  respect 
to  the  natural  frequencies  of  the  individual  cells.  As  before,  har¬ 
monic  forcing  was  applied  throughout  the  frequency  range  of 
1500—3000  Hz,  and  prestress  values  of  150  kN  and  470  1<N  were 
applied. 

Fig.  14  shows  the  maximum  forced  response  amplitude  for  each 
cell  in  the  battery  pack  for  a  prestress  level  of  150  kN,  both  with  and 
without  structural  variations.  All  of  the  amplitudes  in  this  figure  are 
normalized  by  the  maximum  amplitude  of  the  nominal  system.  So 
a  normalized  amplitude  of  1.0  (shown  by  the  blue  horizontal  line) 
represents  the  maximum  response  level  in  the  absence  of  cell-to- 
cell  variations.  For  the  system  with  cell-to-cell  variations 
included,  the  maximum  response  was  first  obtained  for  each  cell 
across  the  entire  frequency  range  of  interest  for  each  Monte  Carlo 
realization.  Then,  the  maximum  value  among  all  10,000  re¬ 
alizations  was  taken.  The  red  bars  plotted  in  Fig.  14  therefore 
represent  the  worst-case  forced  response  amplitudes  for  each  cell, 
and  the  horizontal  red  line  illustrates  the  maximum  response  level 
with  structural  variations.  It  can  be  seen  that  effects  of  cell-to-cell 
variations  lead  to  higher  response  levels  for  all  cells.  The  increase 
is  quite  dramatic  for  some  of  the  cells.  However,  the  increase  for  the 
largest-responding  cell  was  only  about  10%. 

Fig.  15  shows  the  results  for  the  same  simulations  as  in  Fig.  14, 
except  with  the  prestress  level  set  to  470  kN.  Again,  the  struc¬ 
tural  variations  led  to  a  higher  response  for  every  cell.  Further¬ 
more,  there  was  an  increase  of  over  40%  for  the  largest- 
responding  cell.  Therefore,  for  this  prestress  level,  neglecting 
the  cell-to-cell  variations  might  lead  to  a  serious  underprediction 
of  the  maximum  stress,  and  thus  an  overprediction  of  the  battery 
pack  reliability. 
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Fig.  14.  Maximum  normalized  forced  response  amplitudes,  prestress  150  kN. 


As  mentioned  earlier,  when  there  is  no  variation  among  the 
cells,  the  battery  pack  is  a  periodic  structure.  However,  when  the 
cells  have  variation  in  Young’s  modulus,  the  disorder  in  the  struc¬ 
ture  can  lead  to  localization,  in  which  the  vibration  energy  becomes 
concentrated  in  one  or  more  cells.  This  confinement  of  energy  leads 
to  an  increase  in  the  forced  response.  Thus,  it  is  important  to 
investigate  the  effects  of  variations  among  the  cells  in  order  to 
estimate  the  increase  in  maximum  forced  response — as  well  as  the 
attendant  decrease  in  fatigue  life— for  an  HEV  battery  pack. 

5.  Conclusions  and  discussion 

There  is  currently  a  lot  of  research  interest  in  hybrid  electric 
vehicle  (HEV)  battery  performance.  However,  relatively  little 
attention  has  been  paid  to  the  structural  dynamics  of  HEV  battery 
packs.  In  this  paper,  the  structural  dynamic  characteristics  of  HEV 
battery  packs  were  explored  using  a  simple  but  representative 
academic  model  of  a  battery  pack  with  20  pouch  cells.  It  was  shown 


that  the  structural  dynamic  characteristics  are  significantly  affected 
by  the  prestress  level  used  for  joining  cells  within  the  pack.  Also,  it 
was  shown  that  the  dynamic  response  of  the  structure  can  be 
sensitive  to  small  structural  variations  among  the  battery  cells, 
because  the  system  is  a  nominally  periodic  structure  that  features 
high  modal  density. 

Thus,  to  predict  the  fatigue  life  for  such  a  system,  statistical 
calculations  should  be  performed  that  account  for  the  cell-to-cell 
structural  variations.  However,  a  structural  finite  element  model 
of  a  full  HEV  battery  could  easily  have  millions  of  degrees  of 
freedom.  The  large  model  size  makes  it  cumbersome  or  infeasible 
to  run  Monte-Carlo-type  simulations.  Therefore,  parametric 
reduced-order  models  (PROM)  were  derived  to  predict  very  quickly 
the  structural  dynamic  response  of  HEV  batteries  with  random  cell- 
to-cell  structural  variations.  First,  previously  developed  PROM  was 
applied  to  capture  the  different  levels  of  prestress  using  a  numer¬ 
ically  stable  transformation  matrix  and  a  third-order  interpolation 
function.  Second,  to  capture  the  cell-to-cell  variation  in  the  entire 
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Fig.  15.  Maximum  normalized  forced  response  amplitudes,  prestress  470  kN. 
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battery  pack,  new  PROM  was  introduced.  The  new  PROM  is  based 
on  two  key  assumptions:  (1)  mode  shapes  of  the  structure  with 
variations  can  be  represented  as  a  linear  combination  of  mode 
shapes  of  the  structure  with  nominal  parameters,  and  (2)  variability 
in  parameters  in  the  corresponding  cell  can  be  captured  by  mode 
shapes  of  the  nominal  cell  with  its  boundary  displaced  the  same 
amount  as  the  frame. 

As  a  numerical  example,  a  38-DOF  PROM  was  generated  based 
on  the  FEM  of  the  academic  battery  pack  that  has  over  200,000 
DOF.  The  forced  response  results  from  the  PROM  were  found  to 
match  very  well  with  those  from  the  much  larger  FEM  while 
reducing  the  simulation  time  by  a  factor  of  10,000.  The  results  also 
showed  that  prestress  and  small  local  variations  in  the  structural 
parameters  of  the  cells  can  induce  very  large  changes  in  the  global 
response. 

The  main  advantage  of  the  PROM  presented  in  this  work  is  that 
it  can  be  used  to  predict  the  dynamic  response  accurately  relative  to 
the  parent  FEM,  but  at  a  much  lower  computational  cost.  Thus,  it  is 
appropriate  to  use  the  full-order  FEM  as  a  benchmark  for  evaluating 
the  performance  of  the  PROM.  One  can  expect  that  PROM  perfor¬ 
mance  should  ultimately  be  validated  by  using  an  industrial  battery 
pack,  not  using  an  academic  battery  pack.  However,  this  is  an  initial 
study.  Because  the  academic  battery  pack  captures  the  basic  to¬ 
pology  and  key  dynamic  characteristics  of  an  industrial  battery 
pack,  it  is  a  suitable  numerical  example  for  this  initial  investigation. 

Although  this  is  an  initial  study  for  an  academic  example  of  a 
battery  pack,  it  serves  to  illustrate  that  large  battery  packs  with 
many  cells  may  be  nominally  periodic  structures  that  are  sensitive 
to  the  effects  of  small,  random  variations.  In  particular,  they  may 
exhibit  Anderson  localization  and  suffer  from  significant  increases 
in  maximum  forced  response  amplitude  and  stress  levels  in  some 
cells.  In  such  cases,  it  is  important  to  be  aware  of  these  issues  and  to 
analyze  the  effects  of  different  prestress  and  variation  levels  as 
early  in  the  design  process  as  possible.  The  PROM  introduced  in  this 
paper  provides  an  efficient  and  useful  analysis  tool  for  this  purpose. 
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